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ABSTRACT 

We explore the implications of a possible cosmic ray (CR) background generated 
during the first supernova explosions that end the brief lives of massive Population III 
stars. We show that such a CR background could have significantly influenced the 
cooling and collapse of primordial gas clouds in minihaloes around redshifts of z ~ 15 — 
20, provided the CR flux was sufficient to yield an ionization rate greater than about 
10~^^ s~^ near the center of the minihalo. The presence of CRs with energies ^ 10^ eV 
would indirectly enhance the molecular cooling in these regions, and we estimate that 
the resulting lower temperatures in these minihaloes would yield a characteristic stellar 
mass as low as ~ 10 Mq. CRs have a less pronounced effect on the cooling and collapse 
of primordial gas clouds inside more massive dark matter haloes with virial masses 
^ 10® Mq at the later stages of cosmological structure formation around 2; ~ 10 — 15. 
In these clouds, even without CR flux the molecular abundance is already sufficient to 
allow cooling to the floor set by the temperature of the cosmic microwave background. 
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1 INTRODUCTION 

How did the cosmic dark ages end, and how was the homo- 
geneous early Universe transformed into the highly complex 
state that we observe today (e.g. Barkana & Loeb 2001; 
Miralda-Escude 2003; Ciardi & Ferrara 2005)? One of the 
key questions is to understand the first stars, the so-called 
Population III (or Pop III) , since they likely played a crucial 
role in driving early cosmic evolution (e.g. Bromm & Lar- 
son 2004; Glover 2005). For instance, the radiation from the 
first stars contributed to the reionization of the intergalac- 
tic medium (IGM) (e.g. Kitayama et al. 2004; Whalen, Abel 
& Norman 2004; Alvarez, Bromm & Shapiro 2006; Johnson, 
Greif & Bromm 2007), leading to the end of the cosmic dark 
ages. After the first stars exploded as supernovae (SNe) they 
spread heavy elements through the IGM, thereby providing 
the initial cosmic metal enrichment (e.g. Madau, Ferrara & 
Rees 2001; Mori, Ferrara & Madau 2002; Bromm, Yoshida & 
Hernquist 2003; Wada & Venkatesan 2003; Norman, O'Shea 
& Paschos 2004; Greif et al. 2007). Furthermore, this IGM 
metallicity was likely cucial in governing the transition from 
early high-mass star formation to the normal-mass star for- 
mation seen today (e.g. Omukai 2000; Bromm et al. 2001; 
Schneider et al. 2002; Bromm & Loeb 2003; Frebel, Johnson 
& Bromm 2007; Smith & Sigurdsson 2007). 
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SN explosions are also thought to be the site of cosmic 
ray (CR) production (e.g. Ginzburg & Syrovatskii 1969), 
and the CR background built up from the first SNe may af- 
fect the cooling and collapse of primordial gas clouds. Earlier 
studies (Shchekinov & Vasiliev 2004; Vasiliev & Shchekinov 
2006) have shown how the presence of CRs at high redshifts 
could have lowered the minimum mass at which primor- 
dial gas could cool and collapse into virialized dark matter 
(DM) haloes (e.g. Haiman, Thoul & Loeb 1996; Tegmark 
et al. 1997). This earlier study assumes the existence of 
ultra-heavy X particles that decay into ultra-high-energy 
CRs which interact with the cosmic microwave background 
(CMB), leading to the production of ionizing photons. In 
this paper, however, we instead examine CR effects which 
occur through direct coUisional ionization of neutral hydro- 
gen. The free electrons created from the ionization act as a 
catalyst for the formation of H2 (McDowell 1961). Because 
molecular hydrogen emits photons through its rovibrational 
transitions, H2 is able to cool primordial gas at temperatures 
below the threshold for atomic hydrogen cooling (^ 10* K). 
Furthermore, CR ionization can indirectly lead to higher 
abundances of HD, which also acts as a cooling agent in low- 
temperature primordial gas (see Johnson & Bromm 2006). 
However, the presence of CRs can lead to ionization heating, 
as well. Whether this direct heating effect is strong enough 
to counter the additional cooling must be determined and 
will depend on the high-redshift CR energy density. Simi- 
lar to our work, RoUinde et al. (2005, 2006) used models 
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of early star formation to estimate the CR energy density 
in the early Universe, constraining it with the observed Li 
abundances in metal-poor Galactic halo stars. 

In this paper, we will investigate the importance of the 
CR feedback on Pop III star formation by modeling its ef- 
fect on the cooling of primordial gas in two cases: collapse 
within minihaloes, and shocks associated with the virial- 
ization of more massive DM haloes during the later stages 
of structure formation. The outline for this paper is as fol- 
lows. In Section 2 wc discuss CR acceleration and propaga- 
tion in the higli-rcdsliift Universe and how these might differ 
from the present-day CEise. Section 3 describes the evolution 
of primordial gas in minihaloes and in virialization shocks 
when accounting for the effects of CRs. For the minihalo 
case we furthermore discuss how the fragmentation scale 
could change for a sufficiently high CR flux. We present our 
conclusions in Section 4. 



2 COSMIC RAYS IN THE HIGH-Z UNIVERSE 

2.1 Population III star formation 

Though CR effects will be examined for a range of star for- 
mation rates, the typical Pop III star formation rate is taken 
to be that found in Bromm & Loeb (2006) at ^ ~ 15, which 
is approximately 'I', ~ 2 x 10"'^ M0 yr~^ Mpc~^ in a co- 
moving volume. This rate was derived using the extended 
Prcss-Schechter formalism (Lacey & Cole 1993) to model the 
abundance and merger history of cold daxk matter (CDM) 
haloes. In a neutral medium, before the redshift of reioniza- 
tion, Bromm & Loeb (2006) assume star formation occurs 
only in haloes that have become massive enough to enable 
atomic line cooling with virial termparatures above approx- 
imately 10* K. Recent work (see Greif & Bromm 2006) sug- 
gests that these more massive haloes were indeed the dom- 
inant site for star formation. Greif & Bromm (2006) argue 
that about 90 percent of the mass involved in metal-free star 
formation initially cooled through atomic line transitions. 

At higher redshifts such as 2; ~ 20, however, the first 
stars are thought to have formed inside of ~ 10^ M© mini- 
haloes through molecular cooling, and this mode of star for- 
mation is much more significant at this time. For the mini- 
halo case at z ~ 20, Yoshida et al. (2003) estimate the rate 
of star formation through H2 cooling to be ^ , ~ 10~^ M0 
yr~^ Mpc^"*. To account for such differences in these de- 
terminations of star formation rates, we examine a range of 
values spanning multiple orders of magnitude. 

The Pop III initial mass function (IMF) currently re- 
mains highly uncertain, so for this study we do not per- 
form our calculations using a specific IMF. For simplicity, 
we instead assume that those Pop III stars whose masses 
lie in the pair instability SN (PISN) range (140-260 M©) 
have an average mass of 200 M© (e.g. Heger et al. 2003). A 
Pop III initial mass function (IMF) that extends over a large 
range of masses would imply that only a fraction of these 
stars were in the PISN range. Thus, here we assume that 
only slightly less than half of Pop III stars are in this mass 
range, leading to a somewhat more conservative value for 
the CR energy density. Our estimate generally corresponds 
to an IMF peaked around 100 M0 . Due to their high mgiss, 
the first stars had very short lifetimes of about 3 Myr (e.g. 



Bond, Carr & Arnett 1984), and we therefore assume instan- 
taneous recycling of the stellar material. As described below, 
this overall picture of Pop III star formation will be used to 
estimate the average CR energy density in the high-redshift 
Universe. 



2.2 Cosmic ray production 

For this paper CRs arc assumed to have been generated in 
the PISNe that may have marked the death of Pop III stars. 
The CRs are accelerated in the SN shock wave through the 
first-order Fermi process by which high-energy particles gain 
a small percentage increase in energy each time they diffuse 
back and forth across a shock wave. This yields a differential 
energy spectrum in terms of CR number density per energy 
(e.g. Longair 1994): 



dncR 

de 



orm ( ^ \ 

nin \ ^min / 



(1) 



where we will use x = —2 for definitcness, a typical value 
given by Fermi acceleration theory (e.g. Bell 1978a). Here, e 
is the CR kinetic energy, emin is the minimum kinetic energy, 
and ncR is the CR number density. Note that our results are 
somewhat sensitive to the choice of x. Choosing values for x 
that are closer to what is observed in the Milky Way, such as 
X = —2.5 or a; = —3, or using a CR spectrum similar to that 
given in RoUinde et al. (2005, 2006) will increase CR heat- 
ing and ionization, if the overall CR energy density is held 
constant. For these steeper power laws a higher fraction of 
the total CR energy resides in the lower-energy CRs, which 
are the ones that contribute the most to these CR effects, 
as discussed in Sections 3.1 and 3.4 and in the discussion for 
Figure 3. . Such steeper spectral slopes can in fact yield CR 
heating and ionization rates up to an order of magnitude 
higher than for x = —2. Using x = —2 is therefore a more 
conservative choice that does not assume any modifications 
to standard Fermi acceleration theory. 

By equating the total CR energy density Ucr with the 
integral of the differential CR spectrum over all energies, the 
normalizing density factor rinorm is estimated to be 



UCR 



1 ^CR 
10 £min 



(2) 



where we get approximately 1/10 for the coefficient choosing 
c-max = 10^^ eV and emin = 10® eV. The differential energy 
spectrum is therefore 
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with 



This can also be written as 



;7cr(2)«2x10 erg cm (— ) ( in52„,„ 1 



/piSN 



V 0.1 y V 10^2 erg 



(3) 



(4) 



2 X 10-3 MqV \2 X lO-2M0yr-iMpc-3y ' 

where pcR is the fraction of SN explosion energy, JSsn, that 
goes into CR energy, and /pisn is the number of PISNe that 
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occur for every solar mass unit of star- forming material. Here 
we take the above star-formation rate to be constant over a 
Hubble time tn, where 



t«(z)=.2xl0«yr(i±l 



-3/2 



(6) 



evaluated at the relevant redshift. We assume that each star 
quickly dies as a PISN with Ssn = 10®^ erg, appropriate 
for a 200 Mq star (Heger & Woosley 2002), ten percent of 
which is transformed into CR energy (e.g. Rudorman 1974). 
The value of ten percent is derived from Milky Way (MW) 
energetics, and here we have simply extrapolated this to 
PISNe. Very little is known about what value of pcR applies 
to PISNc specifically, so assuming their shock structure to 
be similar to local SNe appears to be a reasonable first guess. 
We choose 1/500 M© for /pisn, so that there is one PISN 
for every 500 Mq of star-forming material. This implies that 
somewhat less than half of the star forming mass falls within 
the PISN range. 

Compared to PISNe, the usual core-collapse SNe (CC- 
SNe) thought to accelerate CRs in the Milky Way have an 
explosion energy that is lower by about an order of magni- 
tude. However, the masses of their progenitor stars arc also 
much lower, ranging from ~ 10 — 4OM0. Thus, if we have an 
IMF extending to this lower mass range, there will be ap- 
proximately an order of magnitude more SNe per unit mass 
of star-forming material. However, because _Esn for low-mass 
(~ 10 Mq) CCSNe progenitors will be lower by a similar 
amount, these two effects may cancel out and result in a to- 
tal UcR that is comparable to the PISN case. Furthermore, 
the difference in shock velocities, Wah, between a PISN and 
a CCSN should not be more than a factor of a few. Though 
CCSNo have a lower explosion energy, this energy is used to 
accelerate roughly ten times less ejected mass, Mej , than in 
a PISN. Simple energy conservation, Esn — ^MejUsh, then 
indicates that the shock velocities of these different SNe will 
be similar. This also leads to similar estimates for emin for 
both cases since the minimum CR energy is expected to 
depend on shock velocity. For a given star formation rate 
the results of our study would therefore change little if the 
source of CRs were CCSNe instead of PISNe. 

It should be pointed out that the explosion mechanism 
of the highest-mass CCSNe progenitors, ~ 40 Mq, is still 
somewhat uncertain and may be associated with very high 
explosion energies comparable to that of PISNe, though if 
most of the CCSNe derive from lower-mass progenitors they 
should still have lower explosion energies as argued above. 
Furthermore, we emphasize that the choice of PISNe versus 
CCSNe is not crucial for our study as long as there exists 
some source of sufficient CR production. These sources do 
not necessarily have to be only PISNe. Constraints on the 
number of PISNe in the early Universe, such as provided 
by the problem of overproduction of metals in the IGM by 
PISNe (e.g. Venkatesan & Truran 2003; Tumlinson et al. 
2004; Daigno ot al. 2006), will therefore not affect the rele- 
vance of this study. 

Finally, it is worth noting that the CR energy densities 
used in this paper are well-within upper limits placed by pre- 
vious studies. For instance, one of the constraints that can 
be used to place an upper limit on Ucb. is the ®Li plateau 
observed in metal-poor Galactic halo stars. While a number 
ratio of ^ Li/H~ 10"" is observed (Asplund et al. 2006), the 



much smaller ratio ^ Li/H~ 10 is predicted by big bang 
nucleosynthesis (BBN). In RoUinde et al. (2005 & 2006), the 
overabundance of ®Li compared to predictions from BBN is 
assumed to have come from CRs generated by the first stars. 
In RoUinde et al. (2005) this assumption was used to derive 
an upper limit on the high-z CR energy density. Assum- 
ing a model with a single burst of Pop HI star formation, 
they found CR energy densities at the time of these bursts 
to range from approximately 10~^° to 10~^^ erg cm~^ for 
burst times ranging from z — 100 to z = 10. Such limits are 
roughly three orders of magnitude larger than the Ucr val- 
ues used in our investigation, so that we do not violate any 
known constraints on CR production in the early Universe. 



2.3 High-redshift GZK cutoff 

The Greiscn-Zatsepin-Kuzmin (GZK) cutoff is an upper 
limit to the energy of a CR if it is extragalactic in origin and 
thus has traveled through the sea of CMB photons (Greisen 
1966; Zatsepin & Kuzmin 1966). This Umit exists because of 
interactions between CRs and CMB photons which lead to 
photo-pion and pair production, reducing the energy of the 
CR for each interaction. The reactions for pair production 
and photo-pion production, respectively, are (e.g. Berezin- 
sky & Grigoricva 1988) 



p"*" -I- 7 ^ 7r+ -I- n , 
p"*" -I- 7 ^ 7r° -h p"*" . 



(7) 
(8) 
(9) 



Photo-pion production is the most important interac- 
tion for CR protons (De Marco 2005). However, this reaction 

can only take place if the energy of a CMB photon is above 
the energy threshold of et = 140 MeV, but this is possible 
in the rest frame of CRs with sufficiently large Lorentz fac- 
tors (e.g. Longair 1994). The reaction will continue until the 
CR energy falls below the corresponding threshold, which 
is eczK = 5 X 10^^ eV in today's Universe. This cutoff has 
recently been observed by the HiRes experiment (Abbasi et 
al. 2007). 

In the high-redshift Universe, however, the GZK cutoff 
will be somewhat lower, as can be seen as follows: While 
in today's Universe the average energy of a CMB i)lioton is 
ecMB = 2.7fcBTcMB = 6 X lO"** eV, at higher redshifts this 
energy will be larger by a factor of (1 -|- ^). In the CR rest 
frame, the CMB photon energy is 



e'cMB »76 X 10-%V(l-h«) 



(10) 



where 7 is the Lorentz factor of a CR proton. Equating ecMB 
with €t, we find 



7 I 



2 X 10" 
~(^ + z) 



(11) 



Wo can now calculate the CR energy for which the threshold 
for photo-pion production is reached: 



£gzk(^;) = 7mHC 



2 3 X 10^° eV 



(1 + ^) 



(12) 



where rrm is the mass of a proton. A more precise calcu- 
lation, carrying out an integration over the entire Planck 
spectrum and over all angles, yields 
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5 X 10^^ eV 
l + z 



2 X 10^* eV 



(13) 



Thus, at redshifts of 10 or 20, the GZK cutoff is around an 
order of magnitude smaller than in today's Universe, giving 
a robust upper limit to the CR energy. 

The GZK cutoff applies only to those CRs that travel 
large distances through the CMB. A CR can undergo a 
maximum of approximately 10 interactions with CMB pho- 
tons before falling below the GZK limit (e.g. Longair 1994). 
Given an interaction cross section of cr = 2.5 x 10~^* cm^ 
for photo-pion production and a CMB density of n-y = 
5 x 10^ cm~^ (1 + zY, the interaction mean-free-path is 



A = ((771^,) 



lO^'^ cm 



(1 



(14) 



The maximum distance from which ultra-high energy CRs 
(UHECRs) with ecR > 5 x 10^^ eV/(l -|- z) could have orig- 
inated is therefore 



lOA : 



10^^ cm 
[l + z)^ 



(15) 



Thus, at « = 20 a UHECR impinging upon a primor- 
dial gas cloud must have been accelerated in a source no 
more than a proper distance of 3 kpc away. The origin of 
UHECRs remains unknown oven in today's Universe, though 
possible sources range from pulsar winds to active galactic 
nuclei (AGN) and gammar-ray bursts (GRBs) (e.g. de Gou- 
veia Dal Pino & Lazaxian 2001; Waxman 2001; Torres et al. 
2002). Though structures at a = 15 and z = 2Q are not yet 
massive enough to form AGN, pulsars and GRBs may be 
plausible UHECR sources at these redshifts (e.g. Bromm & 
Loeb 2006). 



2.4 Magnetic fields 

Magnetic fields are an important component of CR studies 
due to their effects on both CR acceleration and propagation 
through the Universe. For instance, the maximum energy 
a CR can reach with the Fermi acceleration process has a 
linear dependence on the ambient magnetic field strength if 
the growth rate of the CR's energy is limited by its Larmor 
radius (Lagage & Cesarsky 1983). Under this assumption 
strengths of ^ G are necessary to accelerate a CR to 

10^ eV (Zwcibcl 2003). 

The strength, generation, and dispersal of magnetic 
fields at high redshifts, however, is still highly uncertain. One 
of the most widely held views is that magnetic seed fields 
were created soon after the Big Bang and were later am- 
plified in higher-density structures through dynamo mecha- 
nisms. Seed fields are thought to form through various pro- 
cesses including galaxy-scale outflows in the early Universe 
and the Biermann battery mechanism in regions such as 
shock waves and ionization fronts (e.g. Kronberg et al. 1999; 
Gnedin ct al. 2000). For one illustrative example, Ichiki ct 
al. (2006) propose a mechanism by which seed fields are 
created before recombination through second-order cosmo- 
logical perturbations. They calculate the seed field at the 
redshift of recombination, z^cc ~ 1000, to be 



Bo{ZTec) 



lO-^^'G 



A 



10 kpc 



(16) 




Figure 1. Penetration depth as a function of CR energy. We here 
assume a neutral hydrogen number density of nn = 1 cm~^. Only 
clouds with radii larger than about a few hundred pc will entirely 
attenuate the lowest-energy CRs through ionization losses. This 
also implies that low-energy CRs are necessary to contribute sig- 
nificant ionization and heating, as those are the ones that will 
release most of their energy into gas clouds of size ~ 0.1 — 1 kpc. 



where A is the comoving size of the structure in question 
for scales less than about 10 Mpc. If we take into account 
magnetic flux freezing as structures become dense and the 
Universe expands, we get a seed field at virialization of 



Bo(2vir) 



10"^° G 



A 



10 kpc 



2/3 



(1 -I- Z^irf , (17) 



where {ps/ po)z,-„ ~ 200 is the overdensity of the structure in 
question compared to the average density of the Universe at 
Zvir, the redshift at which the structure first virializes. For 
a minihalo with a comoving size of around 10 kpc and viri- 
alization redshift Zvir ~ 20, this gives Bo - 10"^*' G. When 
dynamo effects are taken into account, the magnetic field is 
further amplified and can grow exponentially on a timescale 
determined by differential rotation and turbulence within 
the structure (e.g. Field 1995; Widrow 2002). Magnetic field 
amplification within the accelerating region of a SN remnant 
itself may also increase its strength by up to two orders of 
magnitude (e.g. van der Lann 1962; Bell & Lucek 2001). This 
magnetic field growth can occur through processes such as 
flux freezing in the compressed regions of the SN remnant 
and non-linear amplification through growth and advection 
of Alfven waves generated by the pressure of CRs them- 
selves. Within structures that have already experienced star 
formation, magnetic fields can further be built up through 
field ejection in stellar winds, SN blastwaves, and protostel- 
lar jets (e.g. Machida et al. 2006). 

For these fields to be spread into the general IGM, how- 
ever, there must be a sufficient degree of turbulent mixing 
and diffusion in intcrgalactic regions. While such proccssscs 
arc effective within structures, it is less obvious that they are 
also effective in the IGM at eras soon after star formation 
has begun (Rees 2006). Thus, the build-up and amplificar 
tion of magnetic fields in pristine intcrgalactic matter and 
in the filaments at high redshifts is very uncertain. We can 
estimate the critical magnetic field at z 20 that would 
influence CR propagation by equating the relevant physical 
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distance scale with the Larmor radius tl and finding the 
corresponding magnetic field strength, where for a proton 
rL = jrriHpc^ /{eB), and e is the proton charge. Consider- 
ing the maximum scale a CR could travel at x ~ 20, we write 
rL — f3c/ H{z = 20), where (5c is the CR velocity and H{z) 
the Hubble constant at rcdshift z. For ecR ~ lO'' cV, wc find 
tl — 90 Mpc, corresponding to a magnetic field strength of 
10"^" G, several orders of magnitude higher than what Ichiki 
et al. (2006) estimated for this length scale. When looking at 
CR propagation between minihalocs, the physical distance 
in question is much smaller, tl — 1 kpc. This corresponds 
to a magnetic field strength of 5x10"^'^ G, around an order 
of magnitude stronger than the magnetic field predicted by 
Ichiki et al. (2006) . These seed fields from recombination are 
thus too small to affect CR propagation. 

It is therefore not implausible that at high redshift there 
existed sufficient magnetic fields in the DM haloes to accel- 
erate CRs to relativistic energies while their propagation in 
the IGM and filaments was very close to rectilinear due to 
the weak magnetic fields in these regions. This unconfined 
propagation in the IGM would allow a universal isotropic 
CR background to build up, and this is what we assume in 
our analysis. 



3 EVOLUTION OF PRIMORDIAL CLOUDS 
3.1 Idealized models 

We examine the evolution of primordial gas by adapting the 
one-zone models used in Johnson & Bromm (2006; see also 
Mackey et al. 2003). In the 2 ~ 20 collapsing minihalo case, 
typical halo masses are ~ 10® M©. The minihaloes form 
during hierarchical mergers at velocities too low to cause 
any shock or ionization in the minihalo, so the electron- 
catalysed molecule formation and cooling is insufficient to 
allow stars less than around 100 M© to form (e.g. Bromm, 
Coppi & Larson 1999, 2002; Abel, Bryan & Norman 2002). 
CR ionization and heating, however, can potentially alter 
the thermal and chemical evolution of primordial gas. 

In our model, the minihalo has an initial ionization frac- 
tion of = lO"** and undergoes free-fall collapse. Its den- 
sity therefore evolves according to dn/dt = n/ts, with the 
free-fall time being 



ts ■ 



V 32Gp 



1/2 



(18) 



where p = fimun « mrin and /x = 1.2 is the mean molecular 
weight for neutral primordial gas. The initial density was 

taken to be the density of baryons in DM haloes at the 
point of virialization (e.g. Clarke & Bromm 2003) 



no ^ 0.3 cm 



/i + .y 

V 20 ; 



(19) 



The initial temperature of the gas was taken to be 200 K, 
and the initial abundances were the primordial ones (e.g. 
Bromm et al. 2002). 

As our second case, we consider strong virialization 
shocks that arise in later stages of structure formation dur- 
ing the assembly of the first dwarf galaxies. The correspond- 
ing DM haloes virialize at 2 ~ 10 — 15, and have masses 
ranging from ~ 10* to ~ lO'^" Mq. The conditions during 



the assembly of the first dwarf galaxies show some key differ- 
ences from the former case. The virial velocities of these DM 
haloes axe much greater, implying merger velocites that are 
now high enough to create a shock that can partially ionize 
the primordial gas (Johnson & Bromm 2006 and references 
therein). The post-shock evolution is taken to be roughly 
isobaric (e.g. Shapiro & Kang 1987; Yamada & Nishi 1998). 
We use an initial post-shock temperature of 



Sfca 



(20) 



For an initial density of rips, again found from Equa- 
tion (19), the temperature and density will follow the re- 



lation 



Tn. 



3.2 Thermal and chemical evolution 

In calculating the evolution of the primordial clouds, we 
solve the comprehensive chemical reaction network for all 
the species included in Johnson & Bromm (2006), and con- 
sider cooling due to H, H2, and HD. The temperature of the 
CMB sets the lower limit to which the gas can cool radia- 
tivcly (e.g. Larson 1998). Assuming that CRs with the above 
energy spectrum impinge upon the primordial gas cloud, we 
add the respective heating and ionization rates. Once a low- 
energy CR enters the high-density region of the cloud, each 
time it ionizes an H atom, an electron with average energy 
< >= 35 eV is released (Spitzer & Tomasko 1968). In- 
cluding the ionization energy of 13.6 eV implies that a CR 
proton loses approximately 50 eV of kinetic energy upon 
each scattering. This places a limit on the number of scat- 
terings a CR can undergo in a cloud as well as a limit on 
the distance into the cloud that it can reach before it loses 
all its energy to ionization. This distance can be described 
by a penetration depth 

pee 



where (Schlickeiser 2002) 

- (J) = 1.82 X lO-'eV s-'nHo/(e), 

\Clt / ion 

2/32 



/(e) = (1 -I- 0.0185 ln/3) 



/33 + 2/33' 



and 



1/2 



(21) 

(22) 
(23) 

(24) 



Here, mn is the mass of a proton, /3 = v/c, and (de/dt)ion 
is the rate of CR energy loss due to ionization. The cutoff 
value of Po — 0.01 is appropriate for CRs traveling through 
a medium of atomic hydrogen, since /3oc = 0.01c is the ap- 
proximate orbital velocity of electrons in the ground state of 
atomic hydrogen. When the velocity of CRs falls below Boc, 
the interaction between the CRs and electrons will sharply 
decrease, as will the ionization rate (Schlickeiser 2002), but 
for our study all CRs are assumed to be above this criti- 
cal velocity. Thus, for a given distance D into a cloud, a 
CR has an effective optical depth of D/Dp. Figure 1 shows 
the dependence of Dp on CR energy for a neutral hydrogen 
density of nyo = 1 cm~^, which is typical for densities in 
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Figure 2. Thermal evolution of primordial gas clouds undergoing free-fall collapse inside minihaloes at z = 20 for various combinations 
of Pop III star formation rate and emin- The solid lines show the evolution when the effect of CRs is included, while the dashed lines 
depict the standard evolution where only cooling due to H2 and HD are considered. The dotted line represents the CMB temperature 
floor at 2 = 20. Notice that the presence of CRs can significantly alter the evolution, in particular lowering the temperature, compared 
to the canonical picture of how primordial gas behaves in minihaloes. However, this effect crucially relies on the presence of CRs with 
sufficiently low kinetic energy. 



the minihalo case. Densities will of course greatly increase 
towards the end of the free-fall evolution, but the size of 

the collapsing gas cloud, and thus the distance CRs must 
travel, will decrease. Although the overall attenuation would 
be slightly greater if the time-dependent density evolution 
were accounted for instead of using a constant attenuation 
value, when comparing these two cases the difference is not 
large enough to yield a significant variation in the minihalo's 
temperature evolution. For simplicity, only the typical den- 
sity was therefore used in calculating the attenuation. 

As can be seen in Figure 1, the lowest-energy CRs do 
not get attenuated until they travel a distance of about a 
few hundred pc, so in gas clouds of this size or smaller the 
CR flux will not be significantly attenuated. This also shows 
that the low-energy CRs are the ones that will have the 
greatest ionization and heating contribution to the cloud, 
as they more readily release their energy into the gas. In 
contrast, higher energy CRs will quickly travel through a 
minihalo without transferring much of their energy into the 
gas. They instead lose energy more slowly over much longer 
distances. Accounting for the attenuation yields CR ioniza- 
tion and heating rates of 



and 



50 eV 



(-) 



dncR 
de 



e ' ""de 



(25) 



For 



TIhO Sheat 

These rates can also be written as 



(26) 



and 



CcR(i3) 



c „ in-29 -3 -1 

5 X 10 erg cm s 



2 X 10-is ergcm-3 



5 X 10"^® s'^ 



2 X 10~^5 ergcm~3 



where 



7(e) 



(28) 



(29) 



The factor -Bheat is equal to 6 eV for the minihalo case 
because, though CRs lose about 50 eV of energy after each 
ionization, only about 6 eV of that energy goes toward heat- 
ing in a neutral medium (see Spitzer & Scott 1969; Shull & 
van Steenberg 1985). The value of Shcat increases for me- 
dia with larger ionization fractions due to an increase in 
Coulomb interactions between the newly freed electrons and 
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Figure 3. Minimum temperature reached by tlie minihalo gas versus Ucr (top panel) and spectral slope x (bottom panel). The dotted 
lines denote the CMB temperature, and the dashed line in the upper panel is the upper bound on !7cR based on constraints from the 
^Li measurements. Here emin is kept constant at 10® eV. Note that the cloud cools nearly to the CMB floor for Ucr values ranging over 
almost four orders of magnitude. CR-induecd heating starts to overcome cooling effects and the minimum temperature begins to rise 
only once Ucu values are near the "Li constraint. For shallower spectral slopes, and thus smaller numbers of low-energy CRs for a given 
t^CRi CRs cause little cooling or heating and do not change the minimum temperature of the minihalo gas. 



the medium, and for the ionization fractions typical of the 
virialization shock case, one has E^eat — 26 eV. 

The minimum kinetic energy, emin, with which the CRs 
impinge upon the gas clouds can bo roughly estimated by 
assuming that during Fermi acceleration a particle will ac- 
quire the velocity of the shock wave itself after crossing it 
a single time (sec Bell 1978b). This gives a minimum en- 
ergy €inin ~ imH'ftsh ^ ^V, where Wah ~ 10*kms~^ 
is the PISN shock velocity in its initial blast-wave stage. 
The effects of changing emin will be discussed later. Due 
to the power-law distribution of CR energies, the value of 
the maximum CR kinetic energy, emax, has less importance 
since the number densities and ionization rates of highly rel- 
ativistic CRs are nmch smaller than those of iionrelativistic 
and marginally relativistic ones. We therefore choose a de- 
fault value of emax = 10^® eV, a typical maximum value 
determined from Fermi acceleration theory (e.g. Blandford 
& Eichler 1987). 

Figure 2 shows the gas cloud evolution for the mini- 
halo collapse case for models with and without CR effects 
at z = 20. The heating rate is evaluated at D = 100 pc, the 
characteristic distance to the cloud center, as this is where 
star formation is expected to take place. Various combinar 
tions of = 2 X 10"^ or 2 x 10"^ M© yr"^ Mpc"^ and 
10^ 10^ or 10* oV are considered. The resulting CR 
ionization leads to an increase in the electron abundance 
in the cloud. These electrons act as the only catalyst for 



H2 formation since at this early time in the Universe there 
are no dust grains on which H2 could form. The increased 
electron abundance thus allows for more H2 to form. H2, 
in turn, can also be used in the main reaction that creates 
HD (see, e.g. Johnson & Bromm 2006). With the increase 
in molecular abundance due to CR ionization, the cloud is 
able to cool to temperatures much closer to the CMB floor 
if = 2 X 10"^ M© yr"^ Mpc"^ and emin = 10*^ or lO'^ eV. 
This indirect cooling effect is stronger than the direct CR 
heating effect even for star formation rates up to 100 to 1000 
times higher than those shown here. Futhermore, for 'I', as 
low as 2 X 10"* Mq yr"^ Mpc~^, the CR effects are negligible 
for all emin values explored in this study. 

To further illustrate how the level of CR-induced cool- 
ing can vary with CR parameter values. Figure 3 shows the 
minimum temperature Tmin reached by the minihalo gas ver- 
sus ?7cR and spectral index x given a constant emin = 10® eV. 
It is evident that the cloud can cool to temperatures near the 
CMB floor for f/c;R values ranging over nearly four orders of 
magnitude. At extremely high values of Ucr that are near 
the upper limit set by the ®Li measurements, CR-induced 
heating begins to overcome cooling effects, and the mini- 
halo gas can no longer reach temperatures near the CMB 
floor. Also, for shallower spectral indices, much less of the 
CR energy density is distributed at low energies. For these 
shallower power laws, CRs thus cause only negligible change 
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Figure 4. Thermal evolution of primordial gas clouds experiencing virialization shocks during the assembly of the first dwarf galaxies 
at z = 15. We adopt the same manner of presentation, and the same convention for the lines, as in Figure 2. Notice that the presence 
of CRs has no impact on the evolution, with the exception of extremely high Pop III star formation rates (upper right panel), where 
cooling to the CMB floor is prevented. 



in the minimum gas temperature since they yield almost no 
heating or cooling of the gas. 

Examination of Figure 2 and the corresponding ioniza- 
tion rates for each case show that CRs can facilitate cooling 
in the minihalo to nearly the CMB floor if the ionization rate 
is greater than approximately 10^^'' s~^. Any rate below this 
yields negligible cooling. Considering the dependence of the 
ionization rate on emin and Ucn. in Equations (28) and (29) 
for the minihalo case, we can write 

fc« - ">-" '-' ( .xio-^°:..„-. ) iwk) ■ 

Recalling the dependence of Ucn on redshift and from 
Equation (4), we find that there will be sufficient ionization, 
so that the CR-induced cooling will be significant in the 
minihalo, if the following condition is met: 

(lO-^Moyr-iMpc-s) (^f) (lO^) ^ ^ ' ^^^^ 

Here, we assume that all other variables that determine Ucr 
have the same values used thus far in this study. 

In Figure 4, we show the evolution of a gas cloud in 
the virialization shock case at z = 15. The heating rate is 
now evaluated at D = 500 pc, a typical distance from the 
outer edge of the halo to its densest region. We again exam- 
ine the evolution for various combinations of 4', and emin. 
Note that here we consider star formation rates that are an 
order of magnitude larger than in the minihalo case. Such 



increased rates reflect the later formation times of the first 
dwarf galaxies, so that structure formation has already pro- 
gressed further. In most of these cases the presence of CRs 
has a negligible efi'ect. This is in part due to the enhanced 
CR attenuation resulting from the larger cloud sizes. Fur- 
thermore, unlike in the free-fall minihalo evolution, in the 
case of strong shocks the gas never reaches very high densi- 
ties of neutral hydrogen. Densities for the shocked case are 
instead typically around ^ 0.2 cm~^ early in the halo's evo- 
lution, and 100 cm~^ after the gas temperature has reached 
the CMB floor. The average heating and ionization rates for 
this case are thus lower, but the main explanation for the 
lack of CR-induced difi'erences is that even without CRs the 
shocked region is highly ionized and able to cool through 
molecular transitions to the CMB floor. As shown in John- 
son & Bromm (2006), HD transitions are able to cool the 
gas down to the CMB temperature within a Hubble time. 
Thus, with these cooling mechanisms already in place, CRs 
can only serve to heat the gas, and the CR heating efi'ect 
can dominate over molecular cooling only for more extreme 
star formation rates, as can be seen in the upper right panel 
of Figure 4, where the gas temperature does not reach the 
CMB floor in the presence of strong CR heating. 



3.3 Local CR Feedback 

Our study thus far has assumed that the CRs possibly cre- 
ated in the early Universe all become part of a homogeneous 
and isotropic background. However, if a particular minihalo 
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Figure 5. Thermal evolution of primordial gas in a minihalo under the effect of a PISN explosion 100 pc away (left column) and 10 pc 
away (right column). The PISN is assumed to emit CRs for 2 x 10® yr with a minimum CR energy of 10® eV. We assume that the CR 
flux reaches the cloud during different stages in its collapse, characterized by the corresponding densities: nu = 1 cm~^ (top row) and 
nil = 10^ cm""* (bottom row) . The dashed lines show the evolution without CR flux prsent. The dotted lines show the CMB temperature 
floor at 2: = 20. Notice that the cloud has to be extremely close to the explosion site {d ~ 10 pc) to experience a significant effect. Such 
close proximity, however, is very unlikely due to the strong radiative feedback from the PISN progenitor. 



is within sufficiently close range to a CR-accelerating PISN, 
the flux of CRs from the nearby PISN may have a greater 
effect on the evolution of the minihalo than that from the 
CR background. 

The average CR luminosity associated with the PISN 
is given by 



2 X lO^^rgs" 
\ 1052 erg J 



f PGR ' 

V 0.1 , 



2 X lOByr^ 



-1 



(32) 



where AtsN is the time over which the PISN emits CRs, here 
assumed to be approximately the time that passes from the 
beginning of the PISN to the end of its Sedov- Taylor phase 

of expansion (see Lagage & Cesarsky 1983). We can then 
estimate the CR flux, fen, and energy density, Ucr, emitted 
by the PISN using 

LcR 



fcR 



(33) 



47rd2 ' 

where d is the distance between the CR source and the mini- 
halo. This can also be expressed as 



fcR 



10 erg cm s 



-EsN 



PGR 
0.1 

AtsN 



10^2 erg y \^ 2 X lO^ yr y \^ 100 pc 



The CR energy density can now be estimated to be 
/cR 



UcR'- 



<Pc> 



(35) 



For PGR = 0.1, Es^ = 10^2 erg, and Msn = 2x10-' yr, Ucr 
can then be written as 

/ ^ \-2 r-^^^e'Ae 
t/cR ~ 5 X IQ-i^ergcm-^ — — . (36) 

\100pcy j^"""'' Pe^de 

This energy density is now used in Equations (27) and 
(28) to determine the CR heating and ionization rates due to 
the CR emission from a single nearby PISN. Figure 5 shows 
the thermal evolution of a minihalo under the influence of 
the CR flux from a 10®^ erg PISN at 10 and 100 pc away. 
The PISN was assumed to emit CRs for 2 x 10® yr with a 
minimum CR energy of 10® eV. The resulting CR flux be- 
gins to impinge on the miuilialo at times that correspond 
to two different stages during the collapse, when the gas 
cloud reaches a density of riH = 1 cm~^ and when it reaches 
fiH = 100 cm^"*. At 10 pc, the CR emission allows the gas 
cloud to cool nearly to the CMB floor, while at 100 pc dis- 
tance, it causes only a slight change in the thermal evolution 
of the gas. We also examined the evolution at a distance of 
1 kpc from the PISN, but the CR effects were negligible. 

For such a local burst of CR emission to have a notice- 
able impact on the evolution of the primordial g PISN 
in very close proximity would be required. However, such a 
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Figure 6. Mbe as a function of hydrogen number density for freely-falling clouds inside minihalos at 2 ~ 20. The curves correspond to 
the evolution shown in Fig. 2, and we adopt the same convention for the lines. Notice that the BE mass, which gives a rough indication 
for the final stellar mass, can decrease by up to a factor of 10, given that a sufficient flux of low-energy CRs is present. 



scenario is not very likely. In effect, a cloud at 10 pc distance 
from the explosion would be part of the same minihalo, and 
the strong UV radiation from the PISN progenitor would 
have evaporated any gas in its vicinity by photoionization- 
heating (e.g. Alvarez et al. 2006; Greif et al. 2007). When 
considering realistic cases where the gas evolution in mini- 
haloes could have been significantly impacted by CRs, one 
therefore needs to invoke a universal background that con- 
sists of the global contributions to the CR flux, as opposed to 
the burst-like emission from a single nearby PISN. However, 
at close distances the local CR feedback may still provide 
another source of ionization in nearby gas clouds, indirectly 
leading to increased molecular cooling, thus helping to facil- 
itate collapse and possibly formation of lower-mass stars, as 
discussed below in Section 3.5. 



3.4 Dependence on minimum CR energy 

One of the crucial uncertainties concerning high-redshift 
CRs is the minimum energy £niin of the CRs that impinge 
upon a primordial cloud. Our default value of einin=10^ eV 
derives from a simple estimate for the lowest possible en- 
ergy that a CR proton could gain in a SN shock, though 
other processes may influence this value, possibly increas- 
ing or decreasing it. The minimum CR kinetic energy, how- 
ever, is crucial, because the ionization cross section varies 
roughly as e^p for non-relativistic CRs with kinetic energies 
loss than their rest mass energy but greater than '^lO^ eV, 
the energy corresponding to a velocity of f3o ~ 0.01. Thus, 
higher-energy CRs will travel farther into the cloud before 



first ionizing a particle. Only lower-energy CRs will release 
a large portion of their energy into the cloud, and their ab- 
sence can significantly lower the overall heating and ioniza- 
tion rate inside the cloud. This is especially apparent when 
looking at how the thermal evolution of the minihalo (Figure 
2) changes when emin is increased. For typical stax forma- 
tion rates and Cmin = 10* eV, the impact of CRs becomes 
negligible. It is furthermore interesting to note that for a 
given C/cR and low emin values (~ 10® eV), using a poworlaw 
in momentum instead of a powerlaw in energy will give CR 
heating and ionization rates about an order of magnitude 
smaller than those found in the original calculations. This is 
because using a momentum powerlaw is equivalent to using 
a shallower spectral index at low CR energies, resulting in 
fewer low-energy CRs, and this again illustrates the impor- 
tance of having a sufficient number of low-energy CRs to 
contribute to CR effects. 

Decreasing emin to arbitrarily low values, on the other 
hand, will not yield ever increasing CR ionization rates, as 

the cross section starts to quickly fall off below ^ 10^ oV, 
and our lowest emin of 10* eV is already sufficiently close 
to the peak in ionization cross section to obtain the max- 
imum CR effects. This was confirmed by actually lowering 
the value of Cmin to 1 and 10 keV to see how the heating 
and ionization rates changed. In both the minihalo and viri- 
alization shock case these rates were changed by less than a 
factor of two. 

One of the reasons for the uncertainty in emin is that the 
low-energy part of the CR spectrum is difficult to observe 
even inside the MW. The interaction of CRs with the solar 
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wind and the resulting deflection by solar magnetic fields 
prevents any CRs with energies lower than around 10® eV 
from reaching Earth, and thus MW CRs below this energy 
cannot be directly detected. The Galactic CR spectrum is 
thought to extend to lower energies between 10® and 10^ eV 
(e.g. Webber 1998). However, at non-rclativistic energies the 
CR spectrum in the MW is expected to become much flatter 
due to ionization energy losses that would be much less rele- 
vant in the high- 2; IGM, as even at the redshifts we consider 
the density of the IGM is around three orders of magnitude 
lower than the average MW density. 



3.5 Pragmentation scale 

As is evident in Figure 2, given a strong enough energy den- 
sity, CRs can serve to lower the minimum temperature that a 
collapsing cloud inside a minihalo is able to reach. This could 
have important implications for the fragmentation scale of 
such gas clouds. We estimate the possible change in fragmen- 
tation scale due to CRs by assuming that the immediate pro- 
genitor of a protostar will have a mass approximately given 
by the Bonnor-Ebert (BE) mass (e.g. Johnson & Bromm 
2006) 

Ms.^ 700M,{^f\^^y\ (37) 

where rtt and Tf are the density and temperature of the pri- 
mordial gas at the point when fragmentation occurs. For 
each of the cases shown in Figure 2, the evolution of the 
BE mass was calculated as the cloud collapsed in free-fall. 
This evolution is shown in Figure 6. For the three cases that 
had the most significant CR effects, the respective cooling 
and free-fall times were also examined. When the cooling 
time, tcooi, first becomes shorter than the free-fall time, the 
gas can begin to cool and fall to the center of the DM halo 
unimpeded by gas pressure (Rees & Ostriker 1977; White 
& Rees 1978). As the cloud contracts, however, the density 
increases and the free-fall time becomes shorter, eventually 
falling back below the cooling time. At this point, when 
iff ~ icooi, we evaluate the BE mass. This is the instance 
where the gas undergoes a phase of slow, quasi-hydrostatic 
contraction, what has been termed 'loitering phase' (see 
Bromm & Larson 2004). Slow contraction continues, and 
the gas will leave this loitering regime when its mass is high 
enough to become gravitationally unstable, thus triggering 
runaway collapse. Thus, the 'loitering phase' is the charac- 
teristic point in the evolution when fragmentation occurs 
and the BE mass is relevant, at least to zeroth order, given 
that star formation is too complex to allow reliable predic- 
tions from simple one-zone models such as considered here 
(e.g. Larson 2003). 

Bearing this caveat in mind, we find that, with a suffi- 
ciently strong CR flux present, the gas density at which the 
'loitering phase' occurs is increased by a factor of ^ 10—100, 
and the corresponding temperature decreased by a factor of 
2 to 3. Evaluating Equation (37) shows that the BE mass 
thus decreases by an order of magnitude down to axound 
10 M0. This is the mass scale of what has been termed 'Pop 
11.5' stars (e.g. Mackey et al. 2003; Johnson & Bromm 2006). 
This decrease in fragmentation scale occurs for CR ioniza- 
tion rates greater than ~ 10~^® s"'^. A sufficiently large flux 



of low-energy CRs could thus have facilitated the fragmen- 
tation and collapse of primordial minihalo gas into Pop II. 5 
stars. Again, we emphasize the somewhat speculative na- 
ture of our argument regarding stellar mass scales. However, 
the general trend suggested here, that the presence of CRs 
in the early Universe tends to enable lower mass star for- 
mation, might well survive closer scrutiny with numerical 
simulations. 



4 SUMMARY AND DISCUSSION 

We have investigated the effect of CRs on the thermal and 
chemical evolution of primordial gas clouds in two important 
sites for Pop HI star formation: minihaloes and higher-mass 
haloes undergoing strong virialization shocks. We show that 
the presence of CRs has a negligible impact on the evolution 
in the latter case, since the primordial gas is able to cool to 
the CMB floor even without the help of CR ionization, and 
the direct CR heating is also unimportant unless extremely 
high star formation rates are assumed. Thus, the thermal 
and chemical evolution of haloes corresponding to the first 
dwarf galaxies is rather robust, and the CR emission from 
the deaths of previously formed Pop III stars will not initi- 
ate any feedback in these star formation sites. The impact 
of CRs on gas in minihaloes, which would typically have 
formed before the more massive dwarf systems, could have 
been much more pronounced, given a sufficiently low emin 
and Pop HI star formation rates that are not too low. In 
each case, the effect of direct CR heating is weaker than the 
indirect molecular cooling that follows from the increased 
ionization due to CRs. The additional molecular cooling in- 
duced by the CRs allows the gas in a minihalo to cool to 
lower temperatures. For CR ionization rates above a critical 
value of ~ 10~^^ s~^, we find that the mass scale of metal- 
free stars might be reduced to ^ 10 M©, corresponding to 
what has been termed 'Pop II. 5'. 

If CRs did indeed facilitate Pop II. 5 star formation in 
z = 20 minihaloes, the subsequent production of CRs would 
likely have been signiflcantly reduced. Only one Pop HI star 
is expected to form per minihalo, and it is plausible that 
with CRs present only one Pop II. 5 star would be able to 
form in a given minihalo, as feedback effects of this first 
star may dissipate so much of the minihalo's remaining gas 
that no other stars can form in the same halo. A single 
Pop II. 5 star per minihalo would imply no increase in the 
global number density of stars relative to the case of one 
Pop HI star per minihalo. If Pop II. 5 stars ever formed they 
would all die as lower explosion-energy CCSN. These stars 
would have lifetimes around 10^ yr, not significantly longer 
than for the ~ 100 Mq stars expected to form in minihaloes 
when CRs are unimportant, so Pop II. 5 stars would still 
generate CRs almost instantaneously. However, their con- 
tribution to the CR energy density would be an order of 
magnitude less than that from a PISN due to the reduced 
explosion energy. This reduction in overall CR flux might be 
sufficient to prevent further CR-induced Pop II. 5 star for- 
mation in minihaloes. Stars formed slightly later could then 
again be classical Pop HI stars which would quickly restore 
a high level of CR energy density. The overall impact of CRs 
in the early Universe is thus difficult to determine, owing to 
the intricate feedback between star formation and CR pro- 



© 2007 RAS, MNRAS 000, 1-12 



12 A. Stacy and V. Bromm 



duction. Wo conclude, however, that CRs could significantly 
influence primordial star formation, and it will be important 
to further explore their role with more sophisticated numer- 
ical simulations. 
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